Astronomy & Astrophysics manuscript no. CTA102'vlbi'kin8'mp'LANGED 


©ESQ 2012 


November 16, 2012 





Catching the radio flare in CTA102 
II. VLBI kinematic analysis 

C. M. Fromm 1 , E. Ros 2 - 1 , M. Perucho 2 , T. Savolainen 1 , P. Mimica 2 , M. Kadler 3 , A. P. Lobanov 1 , M. L. Lister 4 , 

Y. Y. Kovalev 51 and J. A. Zensus 1 



(N 

o 

(N 

> 
O 

in 



6 



> 

O 



(N 



X 



Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany e-mail: cfrormn@mpifr.de 
Departament d'Astronomia i Astrofisica, Universitat de Valencia, Dr. Moliner 50, E-46100 Burjassot, Valencia, Spain 
Institut fiir Theoretische Physik und Astrophysik, Universitat Wiirzburg, Am Hubland, 97074 Wiirzburg, Germany 
Department of Physics, Purdue University, 525 Northwestern Avenue, West Lafayette, IN, 47907 USA 
Astro Space Center of Lebedev Physical Institute, Profsoyuznaya 84/32, 117997 Moscow, Russia 



Preprint online version: November 16, 2012 



ABSTRACT 



Context. Very Long Baseline Interferometry (VLBI) observations can resolve the radio structure of active galactic nuclei (AGN) 
and provide estimates of the structural and kinematic characteristics on parsec-scales in their jets. The changes in the kinematics 
of the observed jet features can be used to study the physical conditions in the innermost regions of these sources. We performed 
multifrequency multiepoch Very Long Baseline Array (VLB A) observations of the blazar CTA 102 during its 2006 radio flare, the 
strongest ever reported for this source. These observations provide an excellent opportunity to investigate the evolution of the physical 
properties of blazars, especially during these flaring events 

Aims. We want to study the kinematic changes in the source during the strong radio outburst in April 2006 and test the assumption of 
a shock-shock interaction. This assumption is based on the analysis and modeling of the single-dish observations of CTA 102 (Paper 
I). 

Methods. In this paper we study the kinematics of CTA 102 at several frequencies using VLBI observations. From the modeled 
jet features we derived estimates for the evolution of the physical parameters, such as the particle density and the magnetic field. 
Furthermore, we combined our observations during the 2006 flare with long-term VLBA monitoring of the source at 15 GHz and 
43 GHz. 

Results. We cross-identified seven features throughout our entire multifrequency observations and find evidence of two possible 
recollimation shocks around 0.1 mas (deprojected 18pc at a viewing angle & = 2.6°) and 6.0 mas (deprojected 1 kpc) from the core. 
The 43 GHz observations reveal a feature ejected at epoch f c j = 2005.9 + 0.2, which could be connected to the 2006 April radio flare. 
Furthermore, this feature might be associated with the traveling component involved in the possible shock-shock interaction, which 
gives rise to the observed double peak structure in the single-dish light curves reported in Paper I. 
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1. Introduction 



High-resolution VLBI observations offer the unique possibil- 
ity of studing the s tructural and spectral evolution in AGN on 
parse c-scales (e.g., lLobanov & Zensusl Il999t ISavolainen et al.l 
2002). These observations provide a tool to test several jet mod- 
els and to probe how variations within the first few parsecs con- 
tribute to the high-energy emission. The brightness temperature, 
Th, can be used to derive estimates of the dominant energy- 
loss mechanism and the evolution of intrinsic properties, such 
as the jet speed, the particle density, and the magnetic field (e.g., 
iKadler et al.ll2004l ISchinzel et al.ll2012h . 

The blazar CTA 102 (B2230+114) has been the target of 
several VLBI observations. Within the Radio Reference Frame 
Image Database, the sou rce was observed five times at 8 GHz 
between 1994 and 1998 dPiner et al.ll2007l) . The authors report 
apparent velocities from (-7+ 14) c up to (24+29) c. Higher qual- 
ity images have been obtained within the 15 GHz VLBI obser- 
vations within the VLBA 2 cm-Survey (e.g.. iKellermann et al.l 
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gram (ILister et al 
jet region towards 



2002| and its successor, the MOJAVr£] pro- 
2009a). w an extended 

the southeast spanning a distance up to 
25 mas (deprojecte d 4.5 kpc using a viewing angle of 2.6° 
dJorstad et al.l 120051) ). This region contains several kinks and a 
kinematic analy sis yields apparent speeds between (-3 + l)c 
and (19 + l)c dLister et all l2009bl) . The structure and kine- 
matics of the innermost jet region between March 1998 and 
April 2001 is revealed by the 43 GHz VLBA observations of 
CTA 102 within the Boston University monitoring program^ 
These observations show two stationary features, one close to 
the core (r * 0.1 - 0.2 mas) and the other further downstream 
ate * 2 mas. Furthermore, two newly ejected features (f e j = 
1997.9 + 0.2 and f ej = 1999.5 + 0.1) that are separating from 
the core split into two components further downstream. The flux 
density associated to these moving features outshin e the core for 
a short time in both total and polarized intensity ( Jorsta d~et al.l 



1 Monitoring of Jets in Active galactic nuclei with VLBA 
Experiments http://www.physics.purdue.edu/MOJAVE 

2 http://www.bu.edu/blazars/research.html 
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120051) . The component ejected at t e \ = 1997.9 + 0.2 was con- 
nected to the 1997 flare in CTA 102 dSavolainen et al.ll2002l) . 

iFromm et alJ (1201 ll hereafter Paper I), we analyzed the 
centimeter- and millimeter- wavelength light curves of the his- 
torical flux density outburst of this source that took place in 
2006. Our spectral analysis revealed a double-hump structure 
in the turnover frequency - turnover flux density (v m — S m ) 
plane. The evolution of the physical parameters (e.g., the evo- 
lution of the magnetic field, b (B <x r~ b ), in jet during the 
flare) were extracted from the obtained spectral parameters, 
(v m - S m) by modeling their evolution within the shock-in- je t 
model (Marscher & Gear! Il985t [Bjornsson & Aslaksenl 2000). 
As a possible scenario to explain the observations, we suggested 
an interaction between a traveling and a standing shock wave. 

In this second paper of the series, we test our hypothesis 
of a shock-shock interaction by studying the kinematical evo- 
lution on a parsec scale. Since single-dish observations cannot 
resolve the jet structure, we used several high-resolution multi- 
frequency VLBA observations with a frequency range between 
5 GHz and 86 GHz during the flaring time to study the kinemat- 
ical variations. We combined our observations with the results 
of the 15 GHz MOJAVE observations and additionally analyzed 
the available 43 GHz VLBA observations of CTA 102 within the 
Boston Blazar Monitoring program. This extended data set al- 
lows us to study the morphological evolution of the source over 
nearly two decades (early 1995 until late 201 1). The analysis of 
the core shift and the temporal and spatial spectral evolution will 
be presented in the forthcoming Paper III (Fromm et. al in prep.). 

The organization of the paper is as follows. In Sect. [2] we 
describe the VLBA observations and present the performed data 
analysis. The results of the analysis are presented in Sect. [3]and 
are discussed in Sect [4] A summary and a conclusions are pro- 
vided in Sect. [5] 

Throughout the paper we define the spectral index, a, us- 
ing the relation S v oc y° . The optically thin spectral index, ao, 
can be derived from the spectral slope, s of the relativistic elec- 
tron distribution (N oc E~ s ), via the relation ao = -(s - l)/2. 
We adopt a cosmology with C2 m = 0.27, Qa = 0.73 and Hq 
7 1 kms~' Mpc~'. This results in a linear scale of 8.11 pcmas" 1 
or 26.45 ly mas" 1 for CTA 102 (z=1.037). With these conven- 
tions, 1 mas yr~' corresponds to 52.9 c . 



2. VLBA observations and data analysis 

2.1. Multi-frequency VLBA observations 

We analyzed eight multifrequency VLBA observations with a 
frequency range between 2 GHz and 86 GHz centered around 
the 2006 radio flare. This coverage allowed us to study the flar- 
ing process during its different stages; the rise in the flux density 
from the quiescent stage, the flare, and the decay of the flare. The 
source was observed with all ten antennas of the VLBA, where 
CTA 102 was the main target within the experiment BR 122. In 
the other observations CTA 102 was used as a D-term calibra- 
tor (see Table[T]i. In Table|2]we present the characteristic image 
parameters for the multifrequency VLBA observations. Out of 
the six 86 GHz observations only the May 2005 run has a high 
enough dynamic range to allow us to detect the extended struc- 
ture of the source. Therefore, we cannot provide any kinematic 
results for images at this frequency. 




2004.0 2005.0 2006.0 2007.0 2008.0 2009.0 2010.0 



t[yr] 

Fig. 1. cm-mm single-dish light curves for CTA 102 covering the 
2006 radio flare (see Paper I). The red dashed lines correspond 
to the epochs of multifrequency VLBI observations presented in 
this paper, the black arrows to the 15 GHz, and the gray arrows 
to the 43 GHz VLBA observations. 



Table 1. Used VLBI observations and frequency range 



VLBA-ID 


Date 


Radio bands 11 


Notes 




[yyyy-mm-dd] 






BS157A 


2005-05-19 


CXUKQW 


HN receiver warm 


BS157C 


2005-09-01 


CXUKQW 




BR122A 


2006-04-14 


SCXUKQ 


BR receiver warm 


BR122B 


2006-06-08 


SCXUKQ 


HN 3 hours lost 


BW086B 


2006-10-02 


CXUKQW 


BR focus problems 


BW086C 


2006-12-04 


CXUKQW 




BW086D 


2007-01-26 


CXUKQW 




BW086E 


2007-04-26 


CXUKQW 


SC weather problems 



S=2GHz, C=5GHz, X=8GHz, U=15GHz, K=22GHz,Q=43GHz,W=86GHz 



2.2. Additional 15 GHz and 43 GHz VLBA observations from 
survey programs 

The 15 GHz VLBA observa tions were taken from the MOJAVE 
archive dLister et al.ll2009bh and combined with our 15 GHz ob- 
servations. We analyzed 36 epochs starting from June 1995 until 
October 2010 ((At) = 157 days). 

The 43 GHz VLBA observations were partially taken from 
the Boston blazar monitoring program database 2 . The source 
was on average observed every 58 days between May 2005 and 
January 201 1, which led to a total number of 46 epochs. 
In Figure[TJwe present the temporal correspondence between the 
cm-mm light curves and the VLBA observations of CTA 102. 

2.3. Data reduction 

We used National Radio Astronomy Observatory's (NRAO) 
Astronomical Image Processing System (AIPS) for calibrating 
the data. We performed an amplitude calibration and applied a 
correction to the atmospheric opacity for the high-frequency data 
(v > 15 GHz). The parallactic angle correction was taken into 
account before we calibrated the phases using the pulscal sig- 
nal and a final fringe fit. The time- and frequency-averaged data 

i 1 1 — * I ° 

were imported to DIFMAP (Shepherd 1997), where we used 
the CLEAN algorithm combined with phase and amplitude self- 
calibration and MODELFIT algorithms for imaging and model 
fitting, respectively (see Sect. I2.41 i. 
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Table 2. Image parameters of the multifrequency VLBA obser- 
vations (see Table Q]) 



Frequency 


Epoch 


RMS" 


Speak 


S total 




©min 


P.A. 


[GHz] 


[yyyy-mm-dd] 


[mjy/beam] 


[Jy/beam] 


[Jy] 


[mas] 


[mas] 


[deg] 


5.0 


2005-05-19 


0.21 


1.61 


4.05 


3.65 


1.62 


-8.0 


8.3 


2005-05-19 


0.25 


1.68 


3.91 


2.32 


1.02 


-6.7 


15.3 


2005-05-19 


0.26 


2.22 


3.71 


1.33 


0.55 


-8.8 


22.2 


2005-05-19 


0.55 


2.33 


3.21 


0.95 


0.40 


-15.4 


43.1 


2005-05-19 


0.60 


2.63 


3.16 


0.53 


0.20 


-15.9 


86.2 


2005-05-19 


2.76 


1.94 


2.34 


0.25 


0.11 


-19.3 


5.0 


2005-09-01 


0.32 


1.73 


4.13 


3.36 


1.57 


-3.7 


8.3 


2005-09-01 


0.39 


1.48 


3.23 


1.98 


0.94 


1.4 


15.3 


2005-09-01 


0.31 


1.96 


3.37 


1.21 


0.52 


-6.1 


22.2 


2005-09-01 


0.51 


1.67 


2.49 


0.78 


0.34 


-7.4 


43.1 


2005-09-01 


0.58 


3.48 


4.18 


0.45 


0.18 


-12.3 


2.3 


2006-04-14 


0.51 


2.34 


4.66 


8.60 


3.91 


-7.7 


5.0 


2006-04-14 


0.18 


2.05 


4.48 


3.90 


1.45 


-11.2 


8.4 


2006-04-14 


0.31 


2.94 


4.91 


2.74 


0.97 


-13.2 


15.4 


2006-04-14 


0.26 


4.88 


6.32 


1.46 


0.51 


-14.1 


22.2 


2006-04-14 


0.43 


5.85 


6.84 


1.04 


0.34 


-16.6 


43.2 


2006-04-14 


0.75 


6.43 


7.35 


0.65 


0.15 


-15.8 


2.3 


2006-06-08 


0.29 


2.30 


5.47 


7.75 


3.54 


-2.4 


5.0 


2006-06-08 


0.17 


2.04 


4.40 


3.52 


1.52 


-6.7 


8.4 


2006-06-08 


0.22 


2.78 


4.71 


2.11 


0.95 


-1.1 


15.4 


2006-06-08 


0.25 


4.13 


5.37 


1.26 


0.51 


-8.7 


22.2 


2006-06-08 


0.42 


3.83 


4.76 


0.76 


0.33 


-6.2 


43.2 


2006-06-08 


0.43 


3.62 


4.74 


0.39 


0.18 


-5.0 


5.0 


2006-10-02 


0.22 


2.04 


4.44 


3.51 


1.69 


-9.5 


8.3 


2006-10-02 


0.59 


1.72 


3.45 


2.16 


1.11 


-4.9 


15.3 


2006-10-02 


0.23 


1.96 


2.94 


1.35 


0.60 


-14.7 


22.2 


2006-10-02 


0.46 


1.49 


2.20 


0.76 


0.37 


-7.5 


43.1 


2006-10-02 


0.73 


1.89 


2.78 


0.43 


0.19 


-10.6 


5.0 


2006-12-04 


0.25 


1.84 


4.16 


3.65 


1.75 


-9.7 


8.3 


2006-12-04 


0.23 


1.48 


3.32 


2.11 


1.17 


-4.0 


15.3 


2006-12-04 


0.20 


1.26 


2.23 


1.21 


0.61 


-10.7 


22.2 


2006-12-04 


0.39 


1.08 


1.81 


0.96 


0.55 


9.9 


43.1 


2006-12-04 


0.51 


1.76 


2.40 


0.45 


0.21 


-14.4 


5.0 


2007-01-26 


0.72 


1.43 


3.32 


3.36 


1.71 


-1.6 


8.3 


2007-01-26 


0.26 


1.12 


2.55 


2.06 


1.11 


-0.8 


15.3 


2007-01-26 


0.22 


1.13 


1.99 


1.20 


0.58 


-10.4 


22.2 


2007-01-26 


0.39 


1.22 


1.94 


0.72 


0.39 


-7.7 


43.1 


2007-01-26 


0.66 


2.61 


3.31 


0.38 


0.19 


-9.9 


5.0 


2007-04-26 


0.24 


1.39 


3.59 


3.46 


1.81 


-6.8 


8.3 


2007-04-26 


0.28 


1.16 


2.99 


2.10 


1.14 


-1.2 


15.3 


2007-04-26 


0.20 


1.55 


2.50 


1.23 


0.57 


-9.4 


22.2 


2007-04-26 


0.34 


1.58 


2.24 


0.88 


0.39 


-13.0 


43.1 


2007-04-26 


0.54 


2.54 


3.04 


0.44 


0.20 


-14.7 



" RMS values are determined in a region of the final map without significant emission 



The uncertainties on the obtained fluxes were estimated 
by comparing the total VLBA flux densities with the values 
measured by the University of Michigan Radio Astronomy 
Observatory (UMRAO) program (see Paper I). For the 5 GHz 
and 8 GHz VLBA observations we calculated a flux density dif- 
ference between 10% and 15%, where the difference was be- 
tween 5% and 10% for the 15 GHz ones. The estimates for fre- 
quencies above 15 GHz were more difficult to derive since the 
observations could be easily effected by variations in the wet 
component of the troposphere, but a conservative estimate led 
to an uncertainty valu e of 10% of the UMRAO flux density 
dSavolainen et al.ll2008l) . 



2.4. Model fitting 

To parametrize the jet and trace the evolution of local flux den- 
sity peaks along the jet, we fitted 2D circular Gaussian com- 
ponents to the fully calibrated visibilities using the MODELFIT 
algorithm in DIFMAP. These components were characterized by 
their flux density, S mo d, position r mo d, position angle (PA.), fl mo d 
(measured from north through east) and full width half maxi- 
mum (FWHM), i.e., the deconvolved size of the components. 
We modeled the data at each frequency separately to avoid bias- 
ing effects. Since the number of fitted Gaussians was initially not 
limited, we added only a new component at a given frequency if 
this resulted in a significant decrease in the x 2 value. This ap- 
proach led to a minimum number of Gaussians that can be re- 
ga rded as a rel i able re presentation of the source. As mentioned 
in iLister et all d2009bl) the fitted components do not necessary 
correspond to physical features in the jet and may reflect a math- 
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ematical requirement to properly describe the complex structure 
of the source. To minimize the above-mentioned effect of a fitted 
Gaussian having a purely mathematical nature, we applied the 
following criterion: A fitted feature at a given frequency should 
be traceable at least within five consecutive epochs, and its evo- 
lution in position and flux density should be smooth without any 
strong jumps. If a component satisfies this principle we may in- 
terpret it as a region of locally enhanced emission generated by 
perturbations (shocks or instabilities) in the underlying jet flow, 
the following we identify the most northern component at each 
frequency as the core, e.g., the foot point of the jet and compute 
the position of the additional components relative to the core. 
The uncertainties of the fitt ed components were ty pically around 
5% on the total intensities! Lister & Homan 2005). In this paper 
we follow the suggestion of iLister et al.l ([2009b) that the uncer- 
tainties of the fitted components are around 10% of the compo- 
nent size convolved with the beam size. A more detailed deter- 
mination of the component uncertainties is beyond the scope of 
this paper. 

3. Results 

Here we present the results of modeling the sourc e at differ- 
ent fre quencies. Following the criteria presented bv ILister et al.l 
(2009b), we derive the kinematical parameters, e.g., the apparent 
speed, jSapp, only for those that are identified in five consecutive 
epochs or more, and we regard a component as stationary if the 
angular speed obtained is consistent with /j — within the ob- 
servational accuracy. All analyzed epochs and modeled compo- 
nents at a given frequency are presented in the Appendix. 

3.1. Cross-frequency component identification 

With multifrequency VLBI observations we can study the same 
physical region of a jet at different angular resolutions. In ad- 
dition to the increase in the angular resolution with increas- 
ing frequency, there is the physical effect of the self-absorption 
of the emission. The shape of spectrum of the observed radia- 
tion, assuming that the detected emission is mainly generated 
by an ensemble of relativistic electrons, can be characterized by 
its spectral turnover, e.g., the frequency and flux density at the 
peak (v m and S m ). If the observing frequency is lower than the 
turnover frequency, v,„, the emission will be self-absorbed (opti- 
cally thick), in contrast, if the observing frequency is higher than 
the turnover frequency, the emission will be optically thin. This 
transition between optically thick and optically thin emission oc- 
curs in general in the core region. These two effects are reflected 
in the modeling of the source by splitting the core region into 
several features with increasing frequency. 

For labeling the fitted components we use capital letters for 
the same physical regions and the numbers increase with de- 
creasing core distance. The first criterion for the cross-frequency 
identification of fitted Gaussians is based on their position. Since 
we have not aligned the maps, e.g., corrected for the opacity 
shift, there is a slight shift in the position of these components. 
(The analysis of the opacity shift will be presented in Paper III.) 
However, if a fitted component at a given frequency parametrizes 
a physical region in the jet, the observed angular speed obtained 
for this feature should be consistent with the angular speed de- 
rived from its counterpart at the other frequencies. Based on 
these two criteria we can clearly identify components across the 
different frequencies involved in this study. From our modeling 
we could cross-identify (in frequency and time) four different 
regions, labeled from A to D (see Figs. [2] and [3}. The influence 



of the increased resolution on frequency is clearly visible by the 
amount of detail on the jet structure, which is most obvious on 
the kink around 4 mas away from the core (compare 5 GHz and 
15 GHz image in Fig. |2j. The combined effect of increased res- 
olution and the transition from optically thick and optically thin 
emission can be seen in Fig. [3] 

3.2. Calculation of kinematic parameters 

The motions for the fitted components in CTA 102 can be 
divided into three different categories: i) stationary, ii) non- 
accelearting and iii) accelerating. To derive the kinematics of 
the source, we fit polynomials to the observed x- and y-positions 
of th e fitted components . We follow the procedure presented in 
e.g.. ILister etail (l2009bl) and ISchinzel et all (120121) . using first- 
and second-order polynomials if the components was detected 
in more than ten epochs. (The latter can only be applied to the 
long-term monitoring of the source at 15 GHz and 43 GHz): 



X(t) = fi x (f - f X;0 ) + — (t - f mid ) 



/ \ fly 

y(t) = fly (t - tyfl) + —(t~ ffflid) 



(1) 

(2) 



where f mid = (f min + f max )/2, ? x ,o = fmid - *CW)//*x and f y>0 = 
fmid — y(t m id)/f*y are the component ejection times. Owing to the 
bent structure of CTA 102, the ejection cannot be computed from 
a simple back extrapolation if the trajectory is not clearly ra- 
dial outwards. We only use the results of the second-order poly- 
nomial fit, if there is a significant decrease in the x 2 value as 
compared to the linear fit. Using the average angular speed (fj.) 
obtained, we compute the kinematic parameters, e.g., apparent 
speed, y6 app , and the Doppler factor, 5. The apparent speed of the 
components, f3 app , is derived from the angular speed, fi: 



Aw — 



l+z ' 



(3) 



where Z\ is the luminosity distance and z the redshift. 

The Doppler factor is defined as 6 = r _1 (l - ficos&y 1 , 
where F = (1 - ft 2 )' 1 ' 2 is the bulk Lorentz factor and /3 = v/c. 
Since the /3 and # are not measurable directly from the observa- 
tions, two approaches have been used in the recent past to esti- 
mate the the Doppler factor, 6: 

i) The first one assumes that the observed features are traveling 
at the critical angle, the angle that maximizes the speed of the 
feature. The critical viewing angle, # cr i t can be derived by cal- 
culating the derivative of the apparent speed j6 app = /3 sin - 
/3 cos -&) with respect to & for any given f3. This leads to the re- 
sult that the maximum apparent speed is obtained if cos i9- cl -i t = 

P = [y6a PP /(l +y6app)] ■ For jets seen at this critical angle, the 

Doppler factor is 6 clit = yjl + fil pp . 

ii) We can also use causality arguments and estimate 6 from 
the variab ility of the flux den sity and the size of the component 
(see, e.g.. lJorstad et al.ll2005l) . to get a 'variability' Doppler fac- 
tor 



^var — 



d eS D L 



cAf var (l +z)' 



(4) 



with d e g — 1.8 x FWHM, the effective size of the component 
and Af var = dt/ ln(S 1/52), where S 1 and S2 are the maximum 
and minimum of the component flux density and dt the time dif- 
ference between these two flux density values. The variability 
Doppler factor and the apparent speed can be used to calculate 
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Fig. 2. Uniform weighted VLBA CLEAN images with fitted circular Gaussian components at different frequencies for the 2006 
July observation of CTA 102. The lowest contour is plotted at lOx the off-source rms at 43 GHz and increases in steps of 2. The 
observing frequency and the restoring beam size are plotted above each map. For the labeling we used capital letters for the same 
physical region in the jet and the numbers increase with inverse distance from the core. For a more detailed picture of the core region 
see Fig. [3] 



estimates for the bulk Lorentz factor, F, and the viewing angle, 
■&, according to 



£app + ^var + 1 



26 v 



arctan 



«2 , s2 _ i ' 
Happ ' ^var 1 



(5) 
(6) 



The limited resolution at low frequencies means we cannot re- 
solve some of the components that leads to blending effects be- 
tween two or more features. Under these circumstances the val- 
ues obtained for the apparent speed, y6 app are too low and the 
calculated values for the viewing angle, § are too high. The 
kinematic parameters for the fitted components are presented in 
Sects. EE and [17] Since we observe CTA 102 at 2 GHz only 
for two epochs, April 2006 and June 2006, we neither derived 
kinematic parameters nor present model fitting results in the 
Appendix. 



3.3. Jet region A 

The region labeled as A spreads from r ~ 10 mas to r ~ 18 mas. 
This region could be parametrized with two Gaussian com- 
ponents, labeled as Al and A2, which we traced and cross- 
identified within our 5 GHz, 8 GHz, and 15 GHz observations 
(see Fig. From our multifrequency observations we derived 
an angular speed for Al between (0.15 + 0.09) mas/yr at 5 GHz 
and (0.29 + 0.16) mas/yr at 8 GHz. With the long-term moni- 
toring of CTA 102 at 15 GHz, the angular speed we obtained is 
(0.18 + 0.03) mas/yr. 



The fit of the trajectory of Al by a second-order polynomial 
led to^ d = 5.08, which is no significant increase as compared 
to a first-order polynomial (xL d = 5.14), from which we con- 
clude that there is no detectable acceleration for this component. 
Based on the 15 GHz observations we derived an apparent speed, 
Aipp = (9.4+ 1.9)c (for additional kinematic parameters see Table 
[3}. The angular speed for the second feature in this region, A2, 
spans from (0.11 ± 0.06) mas/yr to (0.21 ± 0.045) mas/yr. We 
could not detect any clear indication of deceleration or accelera- 
tion of this feature from the long-term monitoring at 15 GHz and 
the calculated apparent speeds cover a range between (6 ± 3) c 
and(U±2)c. 

For the presentation of the component trajectories we se- 
lected 15 GHz as the representative frequency. In Fig.|4]we show 
the temporal separation of Al and A2 from the core and the vec- 
tor motion fits are presented in Fig. [5] where the lefthand panel 
shows the position of the components in the sky and the right- 
hand side a zoom of the individual regions with overplotted tra- 
jectory in the contour plot. The components are plotted relative 
to the core position and the contour map corresponds to the 2006 
June observations CTA 102. The lowest contour level is drawn at 
ten times the off-source rms (1 mJy and increase with a factor of 
2). The apparent inward motion of Al (see Fig.|4]i is partly due 
to the nonradial, nearly transversal motion of this feature and in 
contrast, component A2 is mainly moving radially outwards (see 
inlets in Fig. |5j. 

Based on the derived angular speeds, yu, we could, within 
the accuracy and limitation of our time sampling in the multi- 
frequency observations, identifiy Al and A2 across our observa- 
tions. 
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core for all components in region B is plotted in Fig. [6] The vec- 
tor motion of the components Bl, B2, and B3 is presented in 
Fig.0 



15GHz 



2000.0 2005.0 2010.0 

t [year] 

Fig. 4. Temporal separation from the core for the components Al 
and A2 at 15 GHz taken as the representative frequency for the 
evolution of the cross-identified components within this region. 
The solid lines correspond to a first-order polynomial fit to the 
x- and y-position separately. 



3.4. Jet region B 

We parametrized the region further upstream the jet, between 
4 mas and 8 mas with three circular Gaussian components. One 
of these features, labeled as B2, could be identified throughout 
our entire multifrequency data set (5 GHz to 43 GHz). All three 
components, Bl, B2, and B3, show very low pattern speeds. In 
the following we discuss the properties of those components sep- 
arately and in detail. 

The angular speed for Bl is on average comparable to yU = 0. 
However, we found from the 15 GHz observations an angular 
speed of fi = 0.016+0.005 that corresponds to an apparent speed 
y8 ap p = (0.85 +0.25)c. Furthermore, the second-order polynomial 
fit of the y-component of the Bl trajectory leads to a significant 
reduction of the^ 2 value compared to a linear fit, from^ 2 ed = 2.0 
to ^red = 1-0- The value obtained for acceleration /l y = (0.01 ± 
0.002) mas/yr 2 

Similar behavior to Bl is found for B2. The angular 
speed takes values between 0.05 + 0.04 mas/yr and 0.015 + 
0.003 mas/yr, and the corresponding apparent speeds span 
0.77 + 0.14 c to 3 + 2 c. From the long-term monitoring pro- 
grams at 15 GHz and 43 GHz we compute angular velocities 
around 0.015 mas/yr and subluminal apparent speeds y8 app ~ 
0.77 c. Based on the high-resolution 43 GHz observations, we 
could identify B2 within our observational limitations as quasi- 
stationary with a slight radial inward motion (see inlet in Fig.[7j>. 

The innermost component of this triplet is B3. Similar to the 
features farther downstream in region B, it is characterized by 
a small angular speed, where the scatter in // obtained from the 
multifrequency campaign is larger due to the limited time span 
but is, within the uncertainties comparable to the derived angular 
speed from the 15 GHz observation. From these observation we 
calculated an angular speed of fi — (0.027 + 0.006) mas/yr and 
an apparent speed of y8 app = (1.4 + 0.3) c. The trajectory of B3 
shows a nonradial outward motion with slow pattern speed. 

As mentioned before, we used the 15 GHz modelfit results of 
the fitted components in region B as representative of the kine- 
matics within this area. The evolution of the separation from the 




15GHz 



ttyear] 

Fig. 6. Temporal separation from the core for the components 
Bl, B2, and B3 at 15 GHz taken as the representative frequency 
for the evolution of the cross-identified components within this 
region. The solid lines correspond to polynomial fits to the x- 
and y-positions separately. 



3.5. Jet region D 

The third region of the jet spans 1 mas to 3 mas away from the 
core and we parametrized this area with 2 to 3 Gaussians, la- 
beled as Dl, D2, and D3. At these distances from the core, the 
mentioned effects of increasing resolution with frequency lead to 
splitting fitted Gaussians into additional components. Therefore, 
it is difficult to distinguish and clearly cross-identify them. 

The outermost component Dl is located around 3 mas from 
the core and is detected at 8 GHz, 15 GHz, and 22 GHz. We 
found from the long-term monitoring of CTA 102 at 15 GHz in- 
dications of an acceleration in the y direction. This acceleration 
is clearly visible in the y-direction of the 15 GHz representa- 
tive of Dl for t > 2005. Since our time sampling at 8 GHz and 
22 GHz spans 2005.4 < t < 2007.5, the calculated values for the 
angular speed cannot be directly compared to the ones obtained 
at 15 GHz. Therefore, we could not clearly cross-identify this 
feature. The angular and apparent speeds obtained at 15 GHz are 
(0.1 + 0.02) mas/yr and (5.7 + 1.2) c. The component is moving 
on a highly nonradial trajectory outwards, nearly transversal to 
the main jet direction (see inlet in Fig . [9j> ■ 

The feature labeled as D2 is detected throughout our entire 
multifrequency set. As in the case of Dl we found indications 
of an accelerating motion in both, the x and y-component of D2 
from the 15 GHz and 43 GHz observations. As mentioned be- 
fore, owing to the acceleration connected to this component and 
the limited time sampling at 5 GHz, 8 GHz, and 22 GHz the an- 
gular speeds obtained from those observations do not necessar- 
ily correspond to the ones derived from the 15 GHz and 43 GHz, 
but their positions are comparable within the uncertainties. We 
calculated an angular speed /i = (0.016 + 0.001) mas/yr, which 
corresponds to an apparent speed of j8 app = (8.3 + 0.6) c using the 
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long-term monitoring data at 15 GHz and 43 GHz. The trajectory 
of D2 is compared to the main jet direction at to within the un- 
certainties in good agreement with a radial outward motion (see 
inlet in Fig. [9]). 

At 15 GHz we detect one additional feature in this jet re- 
gion, labeled as D3, which we could not distinguish at the other 
frequencies (see Appendix for all fitted components). The com- 
ponent D3 could be interpreted by blending effects in the core 
between two features that could be resolved at higher frequen- 
cies, especially 43 GHz (see following section). 

In Fig. [8] we show the temporal separation from the core for 
Dl, and D2 at 15 GHz and in Fig.|9]we present the vector motion 
fits of these components as crossidentified. 



„ 2.5 
E 



15GHz 



• Dl 

♦ D2 




t [year] 

Fig. 8. Temporal separation from the core for the components D 1 
and D2 at 15 GHz taken as the representative frequency for the 
evolution of the cross-identified components within this region. 
The solid lines correspond to a first-order polynomial fit to the 
x- and y-positions separately. 



3.6. Jet region C 

The cross-identification of the fitted features at distances r < 
1 mas from the core is a difficult task due to both the aforemen- 
tioned splitting of the components and the ejection of new com- 
ponents, which may be associated with flares. These new compo- 
nents can be partially resolved at the highest frequencies but not 
at lower frequencies. In general such a component will appear 
delayed at lower frequencies after it has clearly separated from 
the core and can be resolved at a given frequency. Therefore, 
we concentrate for the core region of CTA 102 on the results 
of the fitting using the 43 GHz observations. We use six circu- 
lar Gaussians, including the core, for the parameterization of the 
core region. Two out of these six components are within our mul- 
tifrequency observations, labeled as CI and C2. 

CI is ejected from the core around 2005.3 with an apparent 
speed of&pp = (11.6 ±3.0) c (ji = (0.22 + 0.06) mas/yr). They- 
component of the CI trajectory shows indication of an nonlinear 
motion but there was not a significant increase in the^ ed to use 
a second-order polynomial fit. Within the uncertainties of our 
observations we classify the motion of CI as radially outwards. 



The second component is ejected from the core in f e j = 
(2005.9 + 0.4) yr with a comparable angular speed to CI ju = 
(0.25 + 0.04 mas/yr corresponds to an apparent speed of yS app = 
(13.0 + 2.1) c. This feature is most probably connected to the 
strong radio outburst in CTA 102 in April 2006 and is moving 
radially outwards (see Sect. |4]i. 

The two additional features, C3 and C4, show a clear non- 
linear trajectory. In the case of C3, the component moves away 
from the core until a distance of 0. 15 mas and keeps this position 
until the last of our observations. A second-order polynomial fit 
reduced the^ ed from 5.0 to 2.3, and the apparent speed obtained 
is/?app = (3.8+0.6) c, which is three times smaller than the appar- 
ent speeds obtained for CI and C2. Owing to the highly curved 
trajectory of the component, we could not compute an estimate 
for the ejection time. The first part of the trajectory could be 
classified as radially outward. However, farther downstream the 
component is moving nonradially inwards. 

The position of C4 remains constant at r ~ 0.1 mas from 
the core until 2009.6 and increases later. The reduced^ 2 value 
decreased from a value of 5.5 to 2.2 by using a second-order 
polynomial instead of a first-order one. This fit leads to an ap- 
parent speed of yS app = (9.1 + 0.5) c. As in the case of C3, we 
could not derive an estimate for the time where the radial sepa- 
ration from the core is consistent with zero (r(f e j) = 0), i. e., the 
ejection time. The trajectory of C3 shows clear indications of a 
nonradial motion. 

In Fig.Q~T]we plot the radial separation from core for all com- 
ponents in the core region and in Fig.[l0]we present the vector 
motion fits. All fitted components in this region have the com- 
mon behavior of having a nearly constant separation from the 
core during their early detection (See inlets in Fig.fTOli. 




2006.0 2007.0 2008.0 2009.0 2010.0 2011.0 2012.0 

ttyear] 

Fig. 11. Temporal separation from the core for the components 
CI, C2, C3, and C4 at 43 GHz taken as the representative fre- 
quency for the evolution of the cross-identified components 
within this region. The solid lines correspond to a first-order 
polynomial fit to the x- and y-positions separately. 



3. 7. Summary of cross-identification 

Here we summarize the results of the previous sections on the 
cross-identification of fitted components. We could clearly iden- 
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tify seven features throughout our entire data set. As mentioned 
before identification of components in the core region is a dif- 
ficult task due to resolution, self-absorption effects and contin- 
uous ejection of new features. However from the 43 GHz ob- 
servations we could detect four new features, were two of them 
where ejected during our multifrequency observations and one 
of them could be connected to the 2006 radio flare (see Sect. @J. 
In Table [3] we present the average kinematic parameters of those 
components. 

Only four components show a clear radial outward motion, 
best seen in CI and C2. These components are detected in the 
core region, where the jet is straight and does not show strong 
bends. On the other hand, the nonradial motion of the majority 
of the cross-identified components may reflect the highly curved 
nature of CTA 102. The inward motion of B2 could be the ef- 
fect of its low pattern speed and the uncertainties in the position 
of the component. In Fig. [12] we show a vector motion map for 
all identified components. In case of the critical angle assump- 
tion the apparent speed corresponds also the lower limit of the 
Doppler factor. 

I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — I 
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15 10 5 

relative RA [mas] 

Fig. 12. Vector motion map for CTA 102. The contour maps cor- 
responds to the stacked 15 GHz observation of CTA 102, and the 
arrows indicate the trajectories of each cross-identified compo- 
nent. The lowest contour level is drawn at 5x the average off- 
source rms (4mJy) and the contours increase in steps of 2. The 
numbers above are the apparent speed of the component. 

The detected acceleration in C3 and C4 could reflect either 
an acceleration of the fluid or the blending effect between a 
stationary feature close to core at r ~ 0.1 mas. However, we 
could not detect indications of any acceleration in CI and C2, 
which are located a t a similar distance from core. Furthermore, 
IJorstad etall (l2005) find in their analysis of earlier 43 GHz VLB! 
observ ations a stationary feature (labeled as Al in Jorsta dTt al.l 
(2005)) at 0.1 mas from the core. If we assume that there is such 
a feature, we could split the components C3 and C4 into a sta- 
tionary and traveling feature and compute the kinematic param- 
eters for the moving ones. Such a procedure leads to an angular 



speed of fx = (0.20 + 0.05)mas/yr 08 app = (10 + 2)c) for C3 
and fx = (0.30 + 0.02)mas/yr 08 app = (15 + l)c) for C4. The 
values obtained are, within the uncertainties, in good agreement 
with the ones for CI and C2, which favors the blending between 
a stationary and traveling components as a possible explanation 
for the trajectories of C3 and C4. 

3.8. Variation in the apparent speed and Doppler factor 

The results of the kinematic analysis of the cross-identified 
components show that the apparent speed increases between 
0.2 mas < r < 1.0 mas (region C in Fig. [3]) from y6 app = 4 to 
/?app = 13. Farther downstream the values decrease from y6 app = 8 
to /?app = 1 between regions D (1.0 mas < r < 5.0 mas) and B 
(5.0 mas < r < 8.0 mas). This fall is related with the very low 
pattern motion of the components in the latter. In the extended jet 
region A (r > 8.0 mas), the apparent speed increases again and 
reaches a constant value (within uncertainties). Figure [T3lshows 
the apparent speed for all cross-identified components versus its 
average position, taken as a representation of its location. The 
apparent speed, /6 app , and the Doppler factor, 6 C iit, assuming the 
components are traveling at the critical angle, § cl -n, are presented 
in Tabled 

The degeneracy of the apparent speed and Doppler factor on 
the intrinsic speed, /3, and the viewing angle, ■&, does not allow 
any direct conclusion on the variation in the fluid speed. This de- 
generacy can be broken if we obtain additional information on 
the Doppler factor from the variability of the flux densities. One 
of the requirements for applying this second method in calculat- 
ing of the Doppler factor is that the flux densities should have 
an ex ponential behavior for several time segments (Jorsta det al.l 
l2005h . The components C2 and Dl fulfill this requirement (see 
Fig. [T4b . and we computed their variability Doppler factor, 6 VSI 
and the corresponding viewing angle, # var . The results of these 
calculations are presented in Table [4] The variability Doppler 
factor decreases from <5 var = 17 ± 3 at (r) ~ 0.4 mas (C2) to 
c\ ar = 8 ± 2 at (r) ~ 2.1 mas (Dl). The viewing angle, # var is 
smaller than the critical angle for both components (see Tables [3] 
and@). We can then conclude that the fluid decelerates between 
(0.4 < r < 2.1) mas. However, there is also a change in the view- 
ing angle, ■&. 



Table 4. Used parameters for the calculation of the variability 
Doppler factor and viewing angle for components C2 and Dl 



component 


d eS [mas] 


Af var 


[yr] 


^var 




C2 (2006.3<t<2007.9) 


0.16 ±0.02 


0.33 ± 


0.06 


17 ±3 


3.2 ±0.2 


Dl(1998.6<t<2003.1) 


0.77 ± 0.05 


3.3 ± 


0.8 


8±2 


6.8 ±0.5 



3.9. Evolution of the brightness temperature T h 
The brightness temperature is defined as 

T h = 1.22 x 10 12 S v /rV 2 (l +z) [K] (7) 

where S v is the flux density in Jy, R the size of the emission 
zone in mas, and v the observing frequency in GHz. Some phys- 
ical conditions, such as the evolution of the magnetic field in 
the source, can be derived from the brig htness temperature (e.g., 
iKadler et all l2004t ISchinzel etal1l2012l) . We present the rela- 
tions between the brightness temperature and the evolution of 
the physical parameters in the framework of the shock-in-jet 
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Comp 



<fi> 
[mas/yr] 



[c] 



^crit "Grit 

in n 



[yr] 



< r > f m ; n f maj[ acceleration 

[mas] [yr] [yr] 



classification 



Al 


0.18±0.()3 


9.4±1.9 


9 


6 




17.0±0.5 


1996.8 


2010.8 


no 


no-radial inward 


A2 


0.17±0.01 


8.7±0.5 


9 


7 


1945±6 


11.2+0.2 


1996.8 


2010.8 


no 


radial outward 


Bl 


0.016±0.005 


0.85±0.25 


1 






7.6±0.3 


1995.6 


2010.8 


yes 


no-radial outward 


B2 


0.015±0.()03 


0.77±0.14 


1 






6.2±0.1 


1996.8 


2010.8 


no 


radial inward 


B3 


0.027±0.006 


1.41±0.32 


2 






5.1±0.2 


1995.9 


2010.8 


no 


non-radial outward 


Dl 


0.10±0.02 


5.7+1.2 


6 


10 




2.1±0.4 


1995.9 


2007.6 


yes 


non-radial outward 


D2 


0.16±0.01 


8.3±0.6 


8 


7 


1997.1±0.4 


1.6±0.4 


2000.0 


2011.0 


yes 


radial outward 


CI 


0.22±0.06 


11.6±3.0 


12 


5 


2005.1 ±0.2 


0.2±0.1 


2005.4 


2006.8 


no 


radial outward 


C2 


0.25±0.04 


13.0+2.1 


13 


4 


2005.9±0.2 


0.4±0.2 


2006.3 


2008.0 


no 


radial outward 


C3 


0.07 ±0.01 


3.8±0.6 


4 


15 




0.2±0.1 


2007.5 


2009.7 


yes 


non-radial outward 


C4 


0.17±0.01 


9.1 ±0.5 


9 


6 




0.2±0.1 


2008.2 


2011.0 


yes 


non-radial outward 
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Fig. 13. Variation in the apparent speed, j6 app along the jet. The 
two red points correspond to the variability Doppler factors, <5 var , 
computed for components C2 and Dl. 



model dMarscher & Gearlll985l) . It is known that the optically 
thin flux density emerging from a relativistic electron distribu- 
tion N = KE S is given by: 

S v ^R 2 x5^ )t2 KB^ l) l 2 V -^l\ (8) 

where R is the width of the jet, x the size of the emitting region, 
B the magnetic field, and s the spectral slope. If the shock prop- 
agates downstream, x is the width of the shock, otherwise x = R. 
The width of the shock front, x, depends on the main energy loss 
mechanis m, namely Compton (1), synchrotron (2), and adiabatic 
losses (3) (Marsch er & Gear|[l985l) and is given by 



Xl cc R-i^K-'B-^'^v- 1 ' 2 (9) 

x 2 cc B -^ 2 6 1/2 v- 1/2 (10) 

x 3 cc R. (11) 

For a more detailed derivation of the equations above, see also 
Paper I. Inserting the equations above into Eq. [8]leads to 

S vA cc R -(s-n»& 6 (s + 4»2 B ( s+ l)/S v - s /2 (12) 

S v , 2 oc R 2 6 (^)/2 KB U-2)/2 y -s/2 (B) 

S v , 3 cc R ^)/2 KB ( S+ W2 v -( S -l),2_ (14) 




1995.0 2000.0 2005.0 2010.0 



t [yr] 

Fig. 14. Variation in the flux densities for the components C2 
(left panel) and Dl (right panel). The gray shaded ares indicate 
the time range used for calculating the variability Doppler factor, 
and the solid line correspond to an exponential fit to the flux 
densities. 



If we assume that the magnetic field, B, the normalization coeffi- 
cient of the relativistic electron distribution, K, and the Doppler 
factor, 5, follow a power law: 

B°cR b KocR k SocR- d , (15) 
Equations [T2TfT4lcan be written in terms of the jet radius 

S^cc/^v 9 ', (16) 
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where the exponent p, includes the dependencies on B, K, and 5, 
and the exponent q-, includes the dependence on the frequency: 



pi = -0s-ll)/8-Z>(,s+l)/8-</(> + 4)/2 (17) 

p 2 = 2 - b(s - 2)/2 - d(s + 4)/2 - k (18) 

Pi = 3-b(s+l)/2-d(s + 3)/2-k (19) 

q\ = s/2 (20) 

qi = s/2 (21) 

q 3 = -(s-l)/2. (22) 



In order to generalize these equations to nonconical jets, the 
distance along the jet, r, is no longer directly proportional to the 
jet radius R. This modification is expressed by R oc r p where 
-1 < p < 1. Finally, we can write the evolution of the brightness 
temperature as 

T b oc r^'v 9 '- 2 , (23) 
where the exponent e, is given by 

e, = p [-(s + 5)/8 - b(s + l)/8 - d(s + 4)/2] (24) 
e 2 = p [-b(s - 2)12 - d(s + 4)/2 - k] (25) 
e 3 = p [1 - 6(s + l)/2 - d(s + 3)/2 - it] . (26) 

The equations above can be simplified further if we assume 
adiabatic expansion k — & a( j = 2(s + 2)/3 and equipartition be- 
tween the magnetic energy density, £ mag oc B 2 and the kinetic 
energy density, £kin ^ X> which leads to b — b eq = (s + 2)/3: 

e ad,eq = ^ r_ (j2 + 5g + n y {6 _ + 4)/2 j (2?) 

^ d - eq = p [-(s 2 +4s + 8)/6 - d(s + 4)/2] (28) 
e* d - eq = p [-(s 2 + 7s + 4)/6 - d(s + 3)/2] . (29) 

These relations can be used to derive the parameter range for 
exponents p, s, k, b, and d using our calculations of the bright- 
ness temperature. Figure [15] shows the variation in 7j along the 
jet for several frequencies, which can be approximated by power 
laws, Tb oc r°, with changes in the exponent in different regions. 

Following the convention of regions defined above, from A 
to D, we derived a decrease in the brightness temperature be- 
tween 0.1 mas and 1 mas (region C). Region D, 1 mas to 4 mas, 
shows a decrease in the Doppler factor (Fig. [TBI and a further 
steepening of the brightness temperature gradient as compared 
to region C. In region B (4 mas to 7 mas) the brightness tem- 
perature values increases and farther downstream (r > 7 mas), 
the brightness temperature decreases again (see Fig.[T5Tl. We ap- 
plied power-law fits to the brightness temperature evolution as 
a function of the core distance to the different regions. The re- 
sults are presented in Tableland are discussed in the following. 
The study of the transversal jet size within the analysis of the jet 
ridge line allowed us to constrain the estimates for the exponent 
p, thereby parametrizing the jet width (Paper III, in preparation). 
The values obtained for the parameter p (jet opening index) are 
presented in Table [6] 

Table 5. Results of the power-law indices from the brightness 
temperature as a function of the distance to the jet basis Tb oc r" . 



region 


C 


D B 


A 


a 


-1.9 + 0.2 


-4.8 ±0.3 2.0 + 0.6 


-2.7 ±0.2 



Table 6. Estimates for the exponent p, defining the jet geometry 
R oc r p (Paper III, in preparation). 



region C D B A 

~~p 0.8 ±0.1 0.8 ±0.1 -1.0 ±0.1 1.3 ±0.2 



We discuss the results in the following terms: i) the jet is in 
equipartition (magnetic energy density is equal to the kinetic en- 
ergy density of the relativistic particles), ii) the magnetic field is 
toroidal (b = 1), and iii) the magnetic field is poloidal (b = 2). 
In the three cases we consider adiabatic expansion and compres- 
sion. We thus assumed that, at distances r > 0.1 mas (depro- 
jected 18pc) Compton and synchrotron los ses can be neglecte d 
as the dominant energy-loss mechanism (Mim ica et al .1 12009). 
The energy losses that the jet undergoes imply a decrease in 
the intensity of the magnetic field and in the Lorentz factors of 
the electrons. This makes it difficult for syn chrotron radiation t o 
dominate the losses far from the jet nozzle (Mi mica et alJ l2009). 
Therefore, we considered only adiabatic losses within adequate 
boundaries for the parameters s, b, k, and d in our analysis. 

3.9.1 . Region C (0.1 mas < r < l.Omas) 

The variation in the brightness temperature in the region B 
(0.1 mas < r < l.Omas) could be best studied at 43 GHz. At 
this frequency we achieved the required resolution to resolve the 
jet, and the long-term monitoring provides us a large number of 
data points that increased the statistical significance. The power- 
law fit gave an exponent of a — —1.9 ± 0.2 (Table|5]l. In Fig. [TBI 
the different panels, labeled 1 to 3, show the range of possible 
parameters for the three different jet configurations proposed (in- 
dicated at the top of each panel). The panels show the variation 
in the exponent that determines the evolution of the Doppler fac- 
tor, 6 oc R d or in terms of distance along the jet 6 oc r~ pd using 
R oc r p in terms of the exponents p and s. We find that for all 
three magnetic field configurations, the Doppler factor should 
increase with distance (d < 0). Among the three cases studied, 
the slowest increase corresponds to the toroidal magnetic field 
(b = 1), and it is faster for the other field configurations, with a 
higher value of abs(d). The increase in the Doppler factor could 
be due to an increase in the fluid speed and/or a decrease in the 
viewing angle. The increase in B app , together with the increas- 
ing 6, suggests an acceleration of the fluid. However, our kine- 
matic analysis of the cross-identified components did not allow 
us to extract any indication of fluid acceleration (see Fig.fTTIand 
Table©. 

3.9.2. Region D (1 mas < r < 4 mas) 

The power law fitted to the evolution of the brightness temper- 
ature {Tb oc r") in this region yields a value of a — -4.8 + 0.3. 
The results of the parameter space study in this case are shown 
in panels 4-6 in Fig. [16] The three different configurations, us- 
ing the value obtained for p (see Table [SJ, lead to two clearly 
different results: On the one hand, if the magnetic field is purely 
toroidal, the Doppler factor had a decreasing exponent through- 
out the entire parameter space. On the other hand, an equipar- 
tition magnetic field and a purely poloidal magnetic filed per- 
mitted a decreasing or increasing Doppler factor, depending on 
the spectral slope. Based on the kinematics obtained, we found 
a slightly decreasing Doppler factor derived for components C2 
and Dl, which is compatible with the toroidal configuration (see 
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panel 5 in Fig.[T6l). but due to the lack of data within this region 
we could not rule out any other possibility. 

3.9.3. Region B (4 mas < r < 7 mas) 

In contrast to the regions C, D, and A, the brightness tempera- 
ture increased with distance in this region. This behavior led to 
a positive value for the exponent a (Tb oc r"). We derived an ex- 
ponent a = 2.0 + 0.6 and at the same distance from the core, we 
found a change from an expanding jet p > to a collimating jet 
p < (see Table [6j. We calculated the possible values for the 
parameter d for a collimating jet p < 0, and the results are pre- 
sented in panels 7-9 in Fig. [16] In the case of a collimating jet, a 
positive value of d corresponds to an increase in the Doppler fac- 
tor, and the opposite is also true (6 oc R~ d &R oc r p — » 5 oc r~ pd ). 
The parameter space study showed that the Doppler factor de- 
creased for all jet configurations within the limits obtained for 
jet expansion rate, p. The major difference between the model is 
that the computed values for d are rather extreme 0.6 < d < 1 .6 
for nontoroidal magnetic field configurations. Assuming that the 
magnetic field is toroidal (fi oc r _I ) in region D, we favor the 
continuation of the magnetic field configurations. 

3.9.4. Region A (r > 7 mas) 

At distances over r > 7 mas, the brightness temperature de- 
creased again with an average exponent of a = -2.7 + 0.2. 
Taking into account that the expansion index of the jet in this 
region is p = 1.3 ± 0.2 (see Tabled, we obtained the possible 
values for the evolution in the Doppler factor (see panels 10 - 
12 in Fig. [TBI , All the studied jet configurations imply an in- 
crease in the Doppler factor and the different models only differ 
in terms of the resulting exponent d. As in the case of region 
C, the slowest increase of the Doppler factors was obtained for 
b — 1 . Considering that the field geometry cannot change from 
toroidal to poloidal, a continuation of a toroidal magnetic field 
configuration is the most plausible one (see panel 1 1 in Fig. [13]). 

4. Discussion 

In the region close to the core (r < 4 mas), we detect sev- 
eral traveling features, most visible at 43 GHz (Figs. [10] and 
fTTb . These features can be interpreted as traveling shock waves, 
which can be g enerated by differe ntial injection pressure at the 
jet nozzle (e.g.. |Perucho et ai1l2008l) . Shock waves travel down- 
stream, compressing the medium ahead of the shock front and 
re-accelerating particles. The increase in radiation induced by 
the shock is observed and can be traced along the jet, until it 
fades due to losses and to the limited dynamical range of the ob- 
servations. In contrast to the expected flux density evolution of a 
traveling shock, we observe that CI, C3, and C4 show some local 
increase in their flux densities while they propagate downstream 
around r ^ 0.1 - 0.2 mas. This could be explained in terms of 
the i nteraction of this shock with a standing r ecollimation shock 
(e.g.. lG6mez et al.lll997UMimica et al.ll2009l) . or by an increase 
in the Doppler boosting of the component if it travels along a 
helic al jet, at the positions where it is closer to the line of sight 
(e.g JHardeefcOOOtlAlov etail 2003). Recollimation shocks may 
appear if the jet and the ambient medium are not in pressure 
equilibrium at the nozzle. Our observations do not allow us to 
resolve any change in the jet direction at this position, and nei- 
ther the Doppler factor (Sect. [3781 and Fig. [TBI nor the brightness 
temperature (Sect. [3791 and Fig. [T5J analysis permit us to obtain 



any clear results from this region, but we did find some indepen- 
dent evidence that favors the presence of a standing shock at this 
position (see Papers I and III, in prep.). 

4.1. The physical parameters along the jet 

4.1 .1 . The Doppler factor 6 

In the previous sections, we have seen that Doppler factor should 
increase in region C (0.1 mas < r < 1 mas), according to 
the shock-in-jet model applied to the observed evolution of the 
brightness temperature. Since our sampling at 43 GHz is not 
dense enough we could not calculate variability Doppler factors 
to test this conclusion. 

The variability Doppler factor, 6 var , computed between re- 
gions C and D shows a decrease with distance (see Fig. [T3land 
Sect. 13.81 ). At the same time, the apparent speed decreases and 
the viewing angle, # var , is smaller than the critical viewing angle, 
# cr i t . This implies that the fluid decelerates within this region. 
This deceleration of the fluid is in good agreement with the cal- 
culated behavior of the Doppler factor based on the gradient in 
the brightness temperature, 7b (see panel 4-6 in Fig.fTol). 

For region B (4 mas < r < 8 mas) we could not derive 
Doppler factors owing to the lack of traveling components. Our 
calculations suggest a decrease in the Doppler factor with re- 
spect to region D, which is independent of the geometry of the 
magnetic field. This deceleration coincides with a region of col- 
limation of the jet (p < 0). Farther downstream, the jet expands 
again and, based on the evolution of the brightness temperature 
within this region, we derived an increase in 5. The obtained evo- 
lution of the Doppler factor can be interpreted within the frame- 
work of an overpressured jet, where the fluid accelerates in the 
regions where the jet expands (p > 0) and decelerates, on aver- 
age, in the regions where the jet reconfines (p < 0). 

4.1.2. The brightness temperature T b 

The brightness temperature (see Sect. |3.9| and Fig.fTSl) decreases 
between 0.1 mas < r < 1 mas (region C), and the decreasing 
trend continues until r - 4 mas (region D). However, there are 
hints of a change in this trend around r - 1 mas, i.e., between 
regions C and D. Between r 4 mas and r 8 mas (region B), 
the brightness temperature increases, and decreases thereafter. 

From the evolution of the brightness temperature in region 
C (0.1 mas < r < 1 mas, see Sect. 13.91 and Fig. [Toll, we de- 
rive that the Doppler factor should be increasing. The value of 
the exponent of 6 as a function of distance depends on the con- 
figuration of the magnetic field, and it increases with b, where 
1 < b < 2. The lowest range for the increase in the Doppler 
factor -0.8 < d < -0.2 is obtained for a toroidal magnetic field 
(b = 1), as seen in panel 2 in Fig.fTol 

In region D ( 1 mas < r < 4 mas) the decrease in the Doppler 
factor could be related to a recollimation and deceleration pro- 
cess, ending up in the stationary component at r - 2 mas. 
dJorstad et al. (2005) report a standing component at the same 
distance.) There is an indication of a slight decrease in the 
growth of the jet radius with distance that could be related to the 
onset of the recollimation process (Paper III, in preparation). The 
highest values for d, implying a steeper decrease in the Doppler 
factor (6 oc R~ d ), are obtained in the case of a toroidal field and 
equipartition. Among these two cases, the latter implies a flatter 
spectral distribution, a situation that will be studied in Paper III. 
The fall in the brightness temperature could be due to the decel- 
eration of the flow in a possible shocked region and to it's not 
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being compensated for by the expected increase in flux density 
at those shocks. 

Between r - 4 mas and r - 8 mas, i.e., region B, the bright- 
ness temperature increases and the jet radius decreases with dis- 
tance (Paper III, in preparation, see Fig. [Tol l. According to our 
calculations, the Doppler factor should decrease. Thus, the in- 
crease in the brightness temperature could also be explained in 
terms of a new process of recollimation of the jet flow, which 
would generate the standing components visible in the region. 
It is difficult to relate the evolution of the Doppler factor with 
distance to the results shown in Fig. [16] owing to the lack of 
data and the large error in the determination of the exponent for 
the evolution of jet radius, p. For larger distances from the core 
the Doppler factor increases, and the brightness temperature in- 
creases again. Within this same interpretation, region A would 
represent a new expansion region. 

4.2. Helical versus overpressured jet 

The viewing angles calculated for the components C2 and Dl 
increase with distance from the core, which indicates a bend in 
the jet away from our line of sight. Since the viewing angles are 
smaller than the critical angles and the apparent speeds decrease, 
we know that the plasma also decelerates. The observed drop in 
the brightness temperature within region D, as compared to re- 
gion C, is a direct consequence of the decreasing Doppler factor 
due to a deviation of the jet away from our line of sight, together 
with the deceleration of the plasma. On the one hand, the sta- 
tionary behavior of components at r ~ 1 .5 mas could indicate the 
presence of a recollimation shock at this position. Such recolli- 
mation shocks are expected to form in overpressured jets, and the 
flow should decelerate when crossing them. On the other hand, 
the increasing viewing angle between regions C and D could in- 
dicate the growth of helical instabilities triggered by pressure 
asymmetries within the jet. 

These instabilities lead to the formation of a helical pattern 
with regions of higher pressure, and hence enhanced emission 
as compared to the rest of the jet flow. Depending on the jet 
properties, we would expect to see the higher pressure region 
tracing helical patterns at higher frequencies and the straighter 
flow at lower frequencies (Perucho et al. 2012). The observed 
morphology of CTA 102 at different frequencies confirms this 
aspect to the extent that the resolution at each frequency per- 
mits. In the case of a strongly over pressured jet, the subsequent 
expansions and collimations would dominate over the linear am- 
plitudes of growing instabilities and we would expect a straight 
flow. The observed changes in the viewing angle imply that the 
jet is not largely over pressured with respect to its surroundings. 
The stability of such a flow and extensions of the regions of en- 
hanced emission depend on the jet properties and the interac- 
tion/ competition between the pinching and helical modes (see, 
e.g JPerucho et al1l2004alfbl) . 

At r > 4 mas from the core, the source is dominated by three 
apparently stationary components at all the observing frequen- 
cies that allow detection. These components also show nearly 
constant size and flux density during the observed period. If the 
components are related to physical regions, this behavior can 
also be explained either by recollimation shocks in a pinching 
jet, by changes in the viewing angle in a helical jet or a combi- 
nation of both. The scenario of an overpressured jet is supported 
by two facts. 1) There are clear hints of jet expansion and subse- 
quent recollimation in the innermost regions (Paper III), which is 
naturally explained in terms of overpressure, and 2) the position 
of stationary components at in region B show that the separation 



between them is growing with distance along the jet. We find 
a separation of ~ 1.0 mas between B3 and B2, and ~ 1.5 mas 
between B2 and Bl. This increase could be naturally caused 
by a decrease in the ambient pressure: If the ambient medium 
is homogeneous, the standing shocks appear at equidistant po- 
sitions and show similar properties, e.g., the jump in pressure 
and density at each shock, whereas if the density in the ambient 
decreases with distance, the shocks appear at increasing separa- 
tions, while pressure jumps become smaller. 

Although the increase in size and the decrease in flux den- 
sity with distance of the stationary components is expected in 
the frame of an expanding jet, this is not an exclusive feature of 
pinching jets, but can also happen in a helical jet. The jet does 
not show any apparent kink between r = 4 mas and r - 12 mas, 
but components Bl, B2, and B3 are slightly misaligned (see 
Fig. which leaves room for short-wavelength kinks in the jet 
in this region, with the changes in the jet direction occurring in 
the short projected distances between components in regions B 
and A. The increase in the separation between components with 
distance to the core could then be explained in this frame if the 
kinks were associated to a helical instability in an expanding jet 
(e.g. JHardee & Hughesll2003tlHardee et aLll200l . The Doppler 
factor seems to grow from region B to region A (see Fig. [T3l . 
However, the lack of components (notice, that Fig. [13] is loga- 
rithmic) does not allow us to discard local changes on this trend, 
i.e., changes in scales of 1 mas. 

It is thus possible to relate, at least partially, the changes in 
the Doppler factor between r — 0. 1 mas and r — 4 mas (regions 
C and D) with changes in the viewing angle, but it is not so easy 
to interpret standing components along the jet within the same 
frame. From the visible changes in the jet direction in the radio 
maps, it seems that the core and the region including compo- 
nents Bl, B2, and B3 would be viewed at short angles than the 
region between D 1 and D2 (2 mas < r < 4 mas) and the region 
beyond B3 (r > 8 mas). These changes could be due to a helical 
pattern in the jet with observed wavelength A \, s 5 mas. The 
misalignment among components B3, B2, and Bl could be ex- 
plained by a small wavelength pattern (A Q ^ 1 mas) developing 
on top of the longer one {A \, s 5 mas). A study of the variations 
in the ridgeline of th e jet with time, such as those performed in 
iPerucho et alj J20l2h . could shed some light on possible helical 
patterns propagating along the jet and on their relation with the 
reported components. 

4.3. The connection to the 2006 radio flare 

As mentioned in Sect. I3.6l we detect several newly ejected fea- 
tures within the core region (see Figs.fTOlandll lb. Furthermore, 
we proposed that the feature at a distance of r ~ 0.1 mas from 
the core is stationary. In Paper I we suggested the interaction be- 
tween a traveling and a stationary recollimation shock as a pos- 
sible process behind the 2006 radio flare and the double hump 
structure in the turnover frequency -turnover flux density plane 
at radio frequencies in CTA 102. In the top panel of Fig.fTTlwe 
show the trajectories of the components CI, C2, C3, and C4. The 
bottom panel shows the flux density evolution of CI, C2, C3, C4, 
and core, together with single-dish light curves at 37 GHz and 
230 GHz, and the evolution of the flux density of the innermost 
components of CTA 102. 

If we assume that there is a standing shock at a distance of 
r ~ 0.1 mas from the core, and if we associate the ejected com- 
ponents with traveling shock waves, there will be shock-shock 
interaction at a given time. This crossing time, t cross , is indicated 
in Fig. \T7\ The derived ejection time for C2, f e j = 2005.9 ± 0.2 
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fits nicely to the first peak in the 230 GHz light curve and to the 
peak in the core flux density (see bottom panel in Fig. [T7b. The 
value obtained for the crossing time between C2 and the pro- 
posed standing shock, ? cr0S s,C2 ~ 2006.3 is in good agreement 
with the peak in the 37 GHz and the second peak in 230 GHz 
single-dish light curve. 

The ejection time for the components C3 and C4 are difficult 
to obtain owing to their highly non-radial trajectory. To calculate 
the ejection times of these components (assuming a blending be- 
tween a traveling and standing shock at r ~ 0. 1 mas), we could 
split the trajectory into moving and standing parts. Since this 
separation could lead to a strong variation in the calculated ejec- 
tion time depending on the position of the separation, we only 
indicate the crossing time for C3 and C4 in Fig[17] However, the 
calculated crossing times f cr0 ss,C3 ~ 2007.7 and f C ross,C4 ~ 2009.3 
correspond to local peaks in the 37 GHz single dish light curves 
(see dashed lines in Fig. IT7T>. 

The peak in the 37 GHz single-dish light curve around 
2008.6 is not connected to the ejection of a new component 
and seems to be connected to C4. However, this region could 
be highly affected by blending and/or misidentification of the 
components. Furthermore, if we assume that there is a stand- 
ing feature at r ~ 0.1 mas, the observed flux density behavior 
could reflect the highly nonlinear interaction of a the shock- 
shock interaction: The strong traveling shock wave associated 
with C2 drags the standing feature downstream. After a given 
time this standing shock wave is re-established at its initial po- 
sition (traveling inward). During this time a new traveling shock 
wave is ejected from the core and collides with the inward trav- 
eling standing shock wave. This proposed scenario could not 
be tested within our observations due to the limited resolution 
within this region and sparse time sampling, and it would require 
high-resolution relativistic hydrodynamic simulations. 

As mentioned before, standing features arise naturally in 
overpressured jets. Therefore, based on the results presented 
above, the overpressured jet model is slightly favored, at least 
in the core region, as a possible jet configuration. 

5. Summary 

In this paper we present the kinematic analysis of the multifre- 
quency VLBA observations of the blazar CTA 102 for eight ob- 
serving epochs during its 2006 radio flare. We combined our data 
with the long-term monitoring at 15 GHz and at 43 GHz pro- 
vided by the MOJAVE and the Boston Blazar Monitoring survey, 
respectively. The modeling of the source at different frequencies 
with several Gaussian features shows that the source consists of 
several apparently stationary components at 4 mas < r < 8 mas 
and a mixture of stationary and traveling ones both in the core 
region and extended regions of the jet. Throughout our entire 
data set, we could clearly cross-identify seven components and 
several newly ejected components in the core region. 

Based on the observed evolution of the cross-identified com- 
ponent trajectories and on the brightness temperature gradients 
obtained, we divided the jet into four different regions. The gra- 
dients in these regions were used to estimate the possible jet con- 
figuration, e.g., orientation of the magnetic field. For our mod- 
eling we calculated the gradients for the evolution of the bright- 
ness temperature according to a generalization of t he shock-in- 
jet m odel, which allows for nonconical geometry (Tiirle r et alj 
2000). We applied three different configurations for the magnetic 
field, b = 1, b = fo eq = + 2)/3, and b = 2 (B oc R~ b ), to the 
observed T/, gradients and derived the evolution of the Doppler 
factor index d (6 oc R~ d ) as a function of the spectral slope, s, and 



the jet geometry. The result of this parameter-space study shows 
that i) in region C there is an increase in the Doppler factor that is 
independent of the magnetic field configuration, ii) the obtained 
decrease in 6 between region C and D favors a jet with a mag- 
netic field toroidal magnetic field and/or in equipartition, iii) the 
decrease in the Doppler factor in this region is due to a deceler- 
ation of the fluid and a change in the viewing angle, iv) for the 
recollimation region B (p < 0) the Doppler factor decreases for 
all magnetic field geometries studied and the smallest decrease 
is obtained for a toroidal field, and v) farther downstream the 
Doppler factor rises again, with the steepnes of the acceleration 
depending on the geometry of the magnetic field. 

The very detailed study of the core region at 43 GHz revealed 
four newly ejected traveling components and one possible sta- 
tionary feature at r ~ O.lmas from the core. We could connect 
one of the components labeled as C2 with the April 2006 radio 
via a linear back-interpolation of its observed trajectory. This 
feature was ejected around 2005.9 from the core with an appar- 
ent speed of (13 + 2) c. The ejection time of C2 fits nicely to the 
first peak in the 230 GHz single dish light curve and the crossing 
time computed for the interaction between C2, and the station- 
ary feature corresponds to the second peak at 230 GHz and to the 
observed global maximum in the 37 GHz single-dish light curve. 

The results were interpreted within the framework of a) an 
overpressured jet, b) a helical jet, and c) a combination of both. 
Within these models the obtained evolution in the Doppler fac- 
tors and brightness temperature are caused by a) an acceleration 
or deceleration of the plasma generated by collimation and ex- 
pansion of the jet, i.e., the pressure mismatch between the jet 
and the ambient medium, by b) a variation in the viewing an- 
gle and the connected Doppler boosting or deboosting and by c) 
nonaxial recollimation shocks that create a transversal pressure 
gradient and lead to a bend in the jet. Based on the analysis of 
the kinematics and the brightness temperature the third option 
was favored. 

To summarize, our results on the component trajectories and 
the brightness temperature evolution suggest a shock-shock in- 
teraction in an overpressured jet as the driving mechanism be- 
hind the 2006 radio flare (see Paper I). However, a detailed anal- 
ysis of the spectral evolution along the jet and, especially, during 
the 2006 radio flare could lead to further evidence of a shock- 
shock interaction. Our hypothesis will then be confirmed or re- 
jected in the last part of the work described in this series of pa- 
pers. 
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Appendix 

Here we show results of modeling for the individual frequencies 
(5 GHz - 43 GHz) during our multifrequency campaign (May 
2005 until April 2007). For each frequency we present contour 
plots of the source, including the position of the fitted circular 
Gaussian components and its structural evolution during our ob- 
servations. Furthermore, we show the temporal evolution of the 
component parameters, e.g., position, flux density and size, and 
the results of a polynomial fit to the cross identified features are 
summarized in Tables 171 -fTTI. The component parameters for all 
frequencies and epochs can be found at the end of the Appendix 
in Tables [l2]and[T3] 



5 GHz VLBI observations of CTA 102 
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Fig. 18. 5 GHz uniformly weighted VLBA CLEAN image 
of CTA 102 observed on 8 June 2006 with fitted circular 
Gaussian components overlaid. The map peak flux density was 
2.0 Jy /beam, where the convolving beam was 3.5 x 1.5 mas at 
P. A. -5.3. The lowest contour is plotted at lOx the off-source 
rms and increases in steps of 2. 
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Fig. 19. Temporal evolution of the separation from the core for 
the 5 GHz components. The color scale corresponds to the flux 
density and the size of the circles to the relative size (FWHM) of 
the components. The solid black lines correspond to a linear fits 
of the component trajectory. 
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Fig. 3. Zoom into the core region for the images presented in Fig. [2] showing the splitting of the C-components with increasing 
frequency. For map details see Fig. [2] 
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Fig. 5. Vector motion fits and position of the fitted component in the sky. Left: 15 GHz CLEAN VLBI map of the 2006 June 
observations of CTA 102. The lowest contour level is drawn at lOx the off-source rms, and the contours increase with steps of 2. 
The observational beam is plotted in the bottom right corner. The blue open symbols indicate the position of the components labeled 
as Al and A2. The filled symbol corresponds to the component position at the middle time, ?o, and the red solid line to the trajectory 
of the feature. The orange squares show the zoom region. Right: Zoom region for the individual components. The dashed gray circle 
are the components size (FWHM) and the dashed black line and black arrow correspond the direction to the core drawn from the 
position of the component at fo. The red solid line illustrates the trajectory of the component, and the direction of the component is 
indicated by the arrow. 



C. M. Fromm et al.: CTA 102 VLBI Kinematic Analysis 




7. Same as Fig.[5]for components Bl, B2, and B3. 
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9. Same as Fig.|5]for components Dl and D2. 
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Fig. 10. Same as Fig.[5]for components CI, C2, C3, and C4, where the core is indicated here by the black pentagon marker 
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Fig. 15. Variation in the brightness temperature, Tf, with distance along the jet for all available epochs and cross-identified com- 
ponents (presented in Tables [T2l and [T3V The nearly vertical variation within 6 mas < r < 8 mas in region B reflect the temporal 
variation in the size of the fitted Gaussians. The solid lines correspond to power-law fits 7b oc r" to the individual regions, and the 
dashed line to a power-law fit to the entire range 0. 1 mas < r < 20 mas, which leads to an exponent a = -1.8 (see text). 
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1 k = K 



= 2 k=K 




Fig. 16. Parameter space for the evolution of the Doppler factor index d as a function of the jet expansion index p and the spectral 

slope s for the different region A - D for three different jet configurations, The first column assumes a adiabatic jet (k = £ ac j = 

2(s + 2)/3) and equipartition between the magnetic energy density and the kinetic energy of the relativistic particles, represented by 
b = b 



eq 



(s + 2)/3. The second and third columns show the evolution of d for an adiabatic jet (k = k a< i) and decreasing magnetic 
field b — 1 or b — 2, respectively. The iso-contours in each panel show the allowed range for the exponent d (6 oc r~ pd ) and the gray 
shaded area the range for the parameter p (R oc r^), with different scales for parameter p in region B. For more details see text. 
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Fig. 17. Evolution of the component separation from the core and the flux density for the innermost components of CTA 102. Top: 
Separation from the core and fitted trajectories for CI, C2, C3, and C4. The gray shaded ares correspond to the ejection epochs for 
CI and C2 and the dashed line to the crossing of a possible stationary feature at r ~ 0.1, mas away from the core (dashed horizontal 
line). Bottom: Flux density evolution of the components including the core and single-dish flux density measurements at 37 GHz 
and 230 GHz. For details see text. 
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Table 7. Results of the kinematic analysis for the fitted components at 5 GHz 



Comp 


H 


/?app 


^min 


^max 


r . 

L mm 




< r > 


^min 


^max 








[mas/yr] 


[c] 


[i] 


[°] 


[1] 


[yr] 


[mas] 


[yr] 


[yr] 


acceleration 


classification 


Al 


0.2±0.1 


8±4 


8 


7 


8 




16.9±0.2 


2005.39 


2007.32 


no 


non-radial outward 


A2 


0.20±0.05 


11±2 


11 


5 


11 


1947+15 


11.6+0.2 


2005.39 


2007.32 


no 


radial outward 


Bl 


0.07±0.07 


4±4 


4 


16 


4 




8±0.3 


2005.39 


2007.32 


no 


non-radial outward 


B2 


0.05±O.O4 


2±2 


2 


21 


2 




6.2±0.1 


2005.39 


2007.32 


no 


non-radial outward 


B3 


0.07±0.06 


4±3 


4 


14 


4 




4.3±0.1 


2005.39 


2007.32 


no 


non-radial outward 


D2 


0.14±0.01 


8±1 


8 


8 


8 


1997.3±0.7 


1.3+0.1 


2005.39 


2007.32 


no 


radial outward 
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Fig. 20. 5 GHz uniform weighted VLBA images of CTA 102 with fitted circular Gaussian components. For better comparison all maps are convolved with a common beam of 
3.5 x 1 .6 mas at P.A. -7.4 and the epoch of the observations is indicated above each contour map. The lowest contour levels is plotted lOx of the maximum off-source rms and 
increases in steps of 2. 
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8 GHz VLBI observations of CTA 102 
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Fig. 21. 8 GHz uniformly weighted VLBA CLEAN image 
of CTA 102 observed on 8 June 2006 with fitted circular 
Gaussian components overlaid. The map peak flux density was 
2.8 Jy /beam, where the convolving beam was 2.1 x 0.9 mas at 
RA. -0.9. The lowest contour is plotted at lOx the off-source 
rms and increases in steps of 2. 
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Fig. 22. Temporal evolution of the separation from the core for 
the 8 GHz modeled features. The color scale corresponds to 
the flux density and the size of the circles to the relative size 
(FWHM) of the components. The solid black lines correspond 
to a linear fit of the component trajectory. 
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Table 8. Results of the kinematic analysis for the fitted components at 8 GHz 



Comp 


H 


/?app 


^min 


^max 






< r > 


^min 


^max 








[mas/yr] 


[c] 


[i] 


[°] 


[i] 


[yr] 


[mas] 


[yr] 


[yr] 


acceleration 


classification 


Al 


0.30±0.20 


15±8 


15 


4 


15 




17.2±0.3 


2005.39 


2007.32 


no 


non-radial outward 


A2 


0.10±0.06 


6±3 


6 


10 


6 




11.5±0.1 


2005.39 


2007.32 


no 


radial outward 


Bl 


0.08±0.08 


4±4 


4 


13 


4 




7.9±0.1 


2005.39 


2007.32 


no 


radial inward 


B2 


0.04±0.03 


2±2 


2 


26 


2 




6.2±0.1 


2005.39 


2007.32 


no 


radial inward 


B3 


0.14±0.07 


7±4 


7 


8 


7 


1970±20 


5.0±0.2 


2005.39 


2007.32 


no 


radial outward 


Dl 


0.8±0.1 


43±7 


43 


1 


43 




2.4±0.5 


2005.39 


2007.32 


no 


non-radial outward 


D2 


0.30±0.02 


16±1 


16 


4 


16 


20()2.3±0.2 


1.2±0.2 


2005.39 


2007.32 


no 


radial outward 
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Fig. 23. 8 GHz uniform weighted VLBA images of CTA 102 with fitted circular Gaussian components. For better comparison all maps are convolved with a common beam of 
2.1 x 1.1 mas at P.A. -2.6 and the epoch of the observations is indicated above each contour map. The lowest contour levels is plotted lOx of the maximum off-source rms and 
increases in steps of 2. 
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15 GHz VLBI observations of CTA 102 
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Fig. 24. 15 GHz uniformly weighted VLB A CLEAN image 
of CTA 102 observed on 8 June 2006 with fitted circular 
Gaussian components overlaid. The map peak flux density was 
4.13 Jy/beam, where the convolving beam was 1.3 x 0.5 mas at 
P.A. -6.2. The lowest contour is plotted at lOx the off-source 
rms and increases in steps of 2. 
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Fig. 25. Temporal evolution of the separation fro m the core for 
the 15 GHz components from lLister et alJd20 09a) and combined 
with our 15 GHz data. The color scale corresponds to the flux 
density and the size of the circles to the relative size (FWHM) 
of the components. The solid black lines correspond to a poly- 
nomial fit of the component trajectory. 
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Table 9. Results of the kinematic analysis for the fitted components at 15 GHz 



Comp 


H 




^min 


$max 


l~min 


'ej 


< r > 


^min 


'max 








[mas/yr] 


[e] 


[i] 


n 


[i] 


[yr] 


[mas] 


[yr] 


[yr] 


acceleration 


classification 


Al 


0.20±0.03 


9±2 


9 


6 


9 




17.0±0.5 


1996.8 


2010.8 


no 


no-radial inward 


A2 


0.20±0.01 


9±1 


9 


7 


9 


1945±6 


11.2±0.2 


1996.8 


2010.8 


no 


radial outward 


Bl 


0.02±0.01 


1.0±0.5 


1 








7.6±0.3 


1995.6 


2010.8 


yes 


no-radial outward 


B2 


0.01±0.01 


0.6±0.5 


1 








6.2±0.1 


1996.8 


2010.8 


no 


radial inward 


B3 


0.03±0.01 


1.4±0.3 


2 








5.1±0.2 


1995.9 


2010.8 


no 


non-radial outward 


Dl 


0.10±0.02 


6±2 


6 


10 


6 




2.1 ±0.4 


1995.9 


2007.6 


no 


non-radial outward 


D2 


0.20±0.01 


8±1 


8 


7 


7 


1997.1±0.4 


1.6±0.4 


2000.0 


2011.0 


yes 


radial outward 


D3 


0.37±0.02 


20±1 


20 


3 


20 




1.1 ±0.4 


2004.3 


2007.6 


no 


non-radial outward 
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Fig. 26. 15 GHz uniform weighted VLBA images of CTA 102 with fitted circular Gaussian components. For better comparison all maps are convolved with a common beam of 
1 .2 X 0.6 mas at P.A. -9.9 and the epoch of the observations is indicated above each contour map. The lowest contour levels is plotted lOx of the maximum off-source rms and 
increases in steps of 2. Notice that the outermost component (Al) for each epoch is plotted as almost overlapping with the jet of the earlier epoch (to its left) in this representation. 
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22 GHz VLBI observations of CTA 102 
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Fig. 27. 22 GHz uniformly weighted VLBA CLEAN image of 
CTA 102 with fitted circular Gaussian components observed 
8th of June 2006 overlaid. The map peak flux density was 
3.8 Jy /beam, where the convolving beam was 0.76 x 0.33 mas 
at RA. -6.2. he lowest contour is plotted at lOx the off-source 
rms and increases in steps of 2. 
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Fig. 28. Temporal evolution of the separation from the core for 
the 22 GHz components. The color scale corresponds to the flux 
density and the size of the circles to the relative size (FWHM) of 
the components. The solid black lines correspond to a linear fit 
of the component trajectory. 
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Table 10. Results of the kinematic analysis for the fitted components at 22 GHz 



Comp 


H 


/?app 


*5min 


$max 


^min 




< r > 


^min 


^max 








[mas/yr] 


[c] 


[1] 


n 


[i] 


[yr] 


[mas] 


[yr] 


[yr] 


acceleration 


classification 


Bl 


0.06±0.11 


3±6 










8.1±0.2 


2005.39 


2007.32 


no 


non-radial outward 


B2 


0.03±0.02 


2±1 


2 


32 


2 




6.2±0.1 


2005.39 


2007.32 


no 


radial inward 


B3 


0.1+0.1 


7±7 


7 


8 


7 




4.9±0.2 


2005.39 


2007.32 


no 


non-radial outward 


Dl 


0.6±0.1 


31±6 


31 


2, 


31 




1.9±0.2 


2005.39 


2006.76 


no 


non-radial outward 


D2 


0.30±0.03 


17±2 


17 
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17 




1.4±0.2 


2005.39 


2007.32 


no 


non-radial outward 
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12 
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Fig. 29. 22 GHz uniform weighted VLBA images of CTA 102 with fitted circular Gaussian components. For better comparison all maps are convolved with a common beam of 
0.8 x 0.4 mas at P.A. -7.6 and the epoch of the observations is indicated above each contour map. The lowest contour levels is plotted lOx of the maximum off-source rms and 
increases in steps of 2. 
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relative RA [mas] 

Fig. 30. 43 GHz uniformly weighted VLB A image of CTA 102 
with fitted circular Gaussian components observed 8th of June 
2006. The map peak flux density was 3.61 Jy/beam, where the 
convolving beam was 0.39 x 0.19 mas at P.A. -5.0. The lowest 
contour is plotted at 5x the off-source rms and increases in steps 
of 2. 



Fig. 31. Temporal evolution of the separation from the core for 
the 43 GHz components. The color scale corresponds to the flux 
density and the size of the circles to the relative size (FWHM) 
of the components. The solid black lines correspond to a poly- 
nomial fit of the component trajectory. 
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Table 11. Results of the kinematic analysis for the fitted components at 43 GHz 



Comp 


H 


/?app 


^min 


^max 


r . 

L mm 


fej 


< r > 


^min 


'max 








[mas/yr] 


[c] 


[i] 


[°] 


[1] 


[yr] 


[mas] 


[yr] 


[yr] 


acceleration 


classification 


CI 


0.20±0.10 


12±3 


12 


5 


12 


2005.1 ±0.2 


0.2±0.1 


2005.4 


2006.8 


no 


radial outward 


C2 


0.25±0.04 


13±2 


13 


4 


13 


2005.9±0.2 


0.4±0.2 


2006.3 


2008.0 


no 


radial outward 


C3 


0.07±0.01 


4±1 


4 


15 


4 




0.2±0.1 


2007.5 


2009.7 


yes 


radial/non-radial outward 


C4 


0.17±0.01 


9±1 


9 


6 


9 




0.2±0.1 


2008.2 


2011.0 


yes 


non-radial outward 




Fig. 32. 43 GHz uniform weighted VLBA images of CTA 102 with fitted circular Gaussian components. For better comparison all maps are convolved with a common beam of 
0.5 X 0.2 mas at P.A. -1 1.9 and the epoch of the observations is indicated above each contour map. The lowest contour levels is plotted lOx of the maximum off-source rms and 
increases in steps of 2. 



Table 12. Fitted components for multifrequency VLBI observations 
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Table 13. Fitted components for multifrequency VLBI observations 
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